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ABSTRACT 

Context. Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step 
in the planet formation process in the core-accretion gas-capture scenario. 

Aims. We aim at providing selected criteria on observational quantities derived from multi-wavelength imaging obser- 
vations that allow to identify dust grain growth and settling. 

Methods. We define a wide-ranged parameter space of discs in various states of their evolution. Using a parametrised 
model set-up and radiative transfer techniques we compute multi-wavelength images of discs at different inclinations. 
Results. Using millimetre and sub-millimetre images we are in the position to constrain the process of dust grain 
growth and sedimentation. However, the degeneracy between parameters prohibit the same achievement using near- to 
mid-infrared images. Using face-on observations in the N and Q Band, the sedimentation height can be constrained. 

Key words, circumstellar matter - planetary systems: proto-planetary discs - radiative transfer - stars: formation 



1. Introduction 

As the existence of extrasolar planets became evident in the 
last century the then established understanding of the for- 
mation of planetary systems also needed serious revision. 
Since then the most accepted evolutionary theory is the 
so-called "core-accretion gas-caption" scenario, featuring a 
build-up of planet cores first and gas-infall on sufficiently 
massive cores later. In the evolution of T Tauri stars, cir- 
cumstellar discs are the natural precursor of a planetary 
system. These discs provide both, dust and gas, from which 
plane ts are expected to be forming (e.g. iKlahr fc Brandnerl 
|2006|) . Key issues in the process of planetary core forma- 
tion are dust grain growth and settling, which start the 
formation process of planetcsimals. 

Dust grains in the interstellar medium (ISM) can be de- 
scribed by a grain size di stribution with a m aximum grain 
size of Qm ax = 250 nm (iMathis et al.l |1977|) . It has been 
shown by Kemper et al.l ( 20041) that amorphous grains, 
which account for more than 98% of the dust in the ISM, 
are smaller than 100 n m after contradicting findings by 
IWitt et all (|200lh and iDraine fc Tanl (|2003h from X-ray 
scattering observations. Molecular clouds are the densest 
regions of the ISM. Circumstellar discs form in the stel- 
lar genesis by conservation of angular momentum from 
a parent interstellar molecular cloud that collapses due 
to instability to its own gravity. Growth of dust grains 
starts already in the molecular cloud a s observed by e.g . 
IWilk ng et all (11980ft and modelled bv lOssenkopfl (119931) 
and Weidenschilling fc Ruzmaikinal (|1994l ). lOrmel etldl 



(|2009j) showed that little coagulation can be expected if the 
cloud lifetime is limited by the typical free fall timescale 
in the order of magnitude of ~ 10 5 yr as the timescale 
of coagulation is of the order ~ 10 7 yr. In contrast, the 



dust grain coagulation times cale in discs is <~ 10 yr (e.g. 
IDullemond fc Dominikll2005t) . 

Evidence for grain growth for T Tauri stars is numerous 
and provided for example by spectro s copy at mid-infrared 
wave lengths (jPrzygodda et alJ 12003: K essler-Silacci et all 
120071 ). Also, low values for the spectral index of the emis- 
sivity in the millimetre regime tracing the densest disc 
regions close to the mid-p lane, /3, indicates grain growth 
(jBeckwith fc Sargentl[l991l ). 

In the evolutionary picture of circumstellar discs, dust 
grains grow from sub-micrometre sizes to millimetre-sizes 
via coagulation. This hit-and-stick process depends on 
numerous properties of the dust grains as for example 
their morphology and their relative velocities. When a 
dust grain has reached a certain size, its coupling to the 
gas will weaken and allow it to settle closer to the mid- 
plane. Effectively, this leads to a m ass-transfer from uppe r 
disc regions towards th e mid- plane (jD'Alessio et al.ll2006[ ). 
IDullemond fc Dominikl (|2004| ) show that as the disc is de- 
pleted in its upper layers, implications to the disc's thermal 
structure and opt ical depth p r ofile a re inevitable. 

According to ICuzzi et al.l (Il993l ). dust grains smaller 
than ~ 1 cm in circumstellar discs can be described within 
the Epstein regime. Small dust grains couple much better 
to their gas environment than large grains. Circumstellar 
discs are understood to be pressure-stabilised against the 
gravity of the star perpendicular to the disc mid-plane. As 
a consequence, larger grains experience a much stronger 
gravitational force along this direction than smaller grains. 
This leads to a settling of large dust grains towards the 
mid-plane. An in-depth treatment of the dust layer in a 
circumstellar disc up t o the onset of hydrod ynamical sta- 
bility can be found in iGaraud fc Lin! p004l ). They argue 
that for a more complete analysis of the situation besides 
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settling, turbulent motion of the gas needs to be considered 
as well. 

Due to the turbulent motion of the gas in deeper disc 
layers dust grains effectively settle in average at maximum 
to their respective sedime ntation height. Compared with a 
disc life time of ~ 10 6 yr (|Sicilia-Aguilar et al.l [20091 ) . sedi- 
mentation is a rather quick process, especially in the upper 
disc layer s. Dust settling depletes these layers in some 10 3 yr 
to 10 4 yr (|Dullemond fc Dominikl I2004D . 

The physical effect that controls the time scale for the 
observables to change is grain growth: As soon as a grain 
has grown to larger sizes, it will settle almost immediately. 

The effects of grain growth on observable quantities of 
circumstellar discs can be seen by examining the luminosity 
of a single dust grain and populations thereof. The lumi- 
nosity of a single grain with radius a, material dependent 
absorption coefficient Qabs and the black-body luminosity 
B(T) of a given temperature T at a certain wavelength A 
can be estimated as 



L x -a 2 Q abSA B\(T). 



(1) 



The total luminosity of a dust population built from these 
grains is then given by 



Q 



abs\ 



J tot x 



(2) 



if one uses the specific material density /3 gr ain, the number 
of grains n, and neglects any temperature dependency of 
Qabs- T he scatter ing properties of dust grains arc affected 
as well (lMidll908l) . 

In addition to the changing optical properties of the 
dust grains grain sedimentation results in the modification 
of the disc density structure. As a consequence, the spatial 
structure of the optical depth t\ changes, yielding different 
temperature profiles in the disc and re-emission features 
of the disc as a function of its evolving state. The altered 
optical depth structure also has implications for the mor- 
phology of scattered light images. 

It is the goal of this study to investigate if the dust 
grain growth and sedimentation induced change of the disc 
structure can be observed. This will provide ideal means to 
verify current dust-coagulation models as well as allow to 
constrain the level of turbulence within circumstellar discs. 
However, the prediction and interpretation of observable 
quantities is rather complex. The optical properties of the 
dust together with the specific density distribution of the 
disc result in a highly non-linear behaviour of the respective 
observables including numerous ambiguities. To cope with 
the complexity of the radiative transfer calculations on the 
model side, Monte Carlo-based codes have proven to be 
quite successful in the last decade. On the observation side, 
the careful interpretation of spectral energy distributions 
(SED) and images taken at different wavelengths offer a 
handle on the nature of circumstellar discs. 

Generally, the SED of more evolved and settled discs 
display a smaller excess at longer wavelengths. This is par- 
allelled in the common classifica tion of classical T Taur i 
stars (CTTS) as introduced by lLada fe Wilkmgl (jl984l) . 
However, the Lada classification evaluates fluxes at wave- 
lengths of 2.2 ism - 10 /j,m while the flux decrease in the 
SED of settled discs occ urs from the near-infrared (N IR) to 
millimetre w a veleng ths (|Dullemond fe Dominikll2005L e.g.). 
M eeus et all (|2009l ) showed that the strength and shape 



of the 10 /Ltm feature in T Tauri stars are related to grain 
growth. They also show that disc flaring and grain sizes 
are slightly related. iD'Alessio et al.l (|2006l ) states essential 
effects of grain growth and dust sedimentation on the SED 
distribution of the disc can be summarised as follows: 

— More settled discs feature a smaller far-infrared (FIR) 
excess. 

— More settled discs yield lower midplane temperatures 
and consequently less flux in the millimetre regime. 

— Discs seen edge-on lack flux at intermediate wavelengths 
as a result of high optical depth at short wavelengths. 
A higher degree of sedimentation yields less extinction 
of the stellar light. 

However, when fitting SEDs to ob servations, model pa 



Thamm et all 11994: 
20091). Further, the 



rameters are often degenerate (e.g. 
iBoss fe Yorkel 119961: ISchegerer etall 
SED is a spatially integrated observable of a disc. Hence, 
it can constrain t hem o nly very weakly as demonstrated by 
iRobitaille et~aTI (|2007l ). 

The effects of dust settling on images of protoplanc- 
tary discs in the optica l wave l ength range are discu ssed by 
iDullemond fc Dominikl (|2004l) . iMcCabe et ail (j2003l) report 
the first resolved image of an edge-on disc in the N band 
and demonstrate that at this wavelength the extent of the 
scattered light nebula is very sensitive to the grain size dis- 
tribution. Experience from modelling individual discs have 
shown that multi- wavelength, spatially resolving images are 
needed to allow constraining of the disc density structure 
and d ust grain prop erties (e.g. Wolf et aTl 20031 : iPinte et all 
120081: ISauter et all 120091: IDuchene et al.ll201(tl . While mil- 
limetre observations are sensitive only to radiation being 
emitted from dust in the densest region within the disc, NIR 
images are dominated by scattered stellar light from dust in 
the circumstellar envelope and the upper, optically thin disc 
layers. These observations trace different physical processes 
(scattering/re-emission) in different regions of the circum- 
stellar environment, but they are both strongly relate d to 
the dust properties in the system. IWatson et al.l (|2007l ) dis- 
cuss the interpretation of images of multi-wavelength ob- 
servations in detail. For example, the disc mass and the 
dust opacity cannot be disentangled as only their product 
is observable. Effects of parameters such as the disc dust 
mass or the inner and outer radii on images have been 
studied, both mo d elling single objec t s and compilations 
(IPinte et al] 120081: ISauter et alj [2001 iBrown et all 120091: 
ISicilia-Aguilar et al.l2007f) . Stratified disc models with large 
dust grains close to the midplane and small grains in the 
upper layers and an envelope have been used to repro 



duce o bservations o f sever a l discs such as IDuchene et al 
(|2010 h : ISauter et al.l (|2009l) : IPinte et all (l2008h : IWolf et al 
20031 ). In these studies a multi- wavelength approach is 



used in order to constrain disc layers with different grain 
sizes. Distinct objects have been modelled with different 
dust models. For instance, some aut hors include grap hite 
(IWolf et al.ll200i and some do not (jPinte et al.ll2008D . In 
the case of CB 26 , this choice has im plications for the max- 
imum gram size (jSauter et al.ll2009D . 

However, predictive studies analysing the impact of 
grain growth and sedimentation on observable quantities in 
a large wavele ngth range have only considered the effects on 
the S EP (e.g. ID'Alessio et al.|[2006HDullemond fc Dominikl 
I2004T) so far. 
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As the next generation of telescopes and interferometers 
with high spatial resolution will commence in the not so 
distant future, it is vital to investigate whether and how 
grain growth and settling can be constrained by spatially 
highly resolved images of circumstellar discs. 

The present paper provides a first approach in address- 
ing the need for a general investigation of the effects of 
dust grain growth and sedimentation on images from the 
near-infrared to the millimetre-regime. We identify indica- 
tors that allow the classification of observations and point 
out directions of interest that are worth a more detailed 
approach in subsequent studies. 

Addressed key questions are: 

1 . What are the general effects of the grain growth induced 
mass relocation on images of a disc seen at different 
inclinations and at different observing wavelengths? 

2. Which tracers of dust grain growth and sedimentation 
are provided by edge-on millimetre images in general? 

3. What can be learned by combining mid-infrared images 
at various inclinations? 

4. How can grain growth and sedimentation be constrained 
by near-infrared images? 

The present study will necessarily invoke several simplifica- 
tions. Modelling of images of an individual object will pro- 
vide a better insight in its particular evolutionary stage. 
Nevertheless, this study aims at providing a first, non- 
comprehensive guide into the field of imaging dust grain 
growth and sedimentation. 

This paper is organised as follows: Section [2] introduces 
the model, the parameters and the methodology employed. 
In Section [3] we illustrate the general effects of dust grain 
growth and sedimentation in images for a fiducial model. 
We also present tracers for grain growth and sedimentation 
in millimetre and mid-infrared images that are distinctive 
from other disc parameters. In Section @] we discuss cases 
for which this fails. 



2. Method 

To investigate the effects of dust grain growth and sedimen- 
tation on observable quantities, four steps are needed: First, 
a model of a circumstellar disc has to be set up that in- 
cludes the effects of dust settling and coagulation. Secondly, 
for all model parameters it has to be determined, whether 
they shall be set to some typical value or if they should 
be allowed to vary within a certain range. In the third 
step, radiative transfer simulations for every model have 
to be performed. As a concluding fourth step, the results of 
the radiative transfer simulations have to be investigated. In 
this Section, we give a description of the first three steps. 

2.1. The model 

The basis of all our models is a parametrised disc density 
distribution. The dust density distribution pdisc of such a 
disc can be written as 



Pdisc(^) = Po 





' I 


' z~ 




UJ exp ( 




X 


1 









(3) 



where z is the usual axial coordinate with z = correspond- 
ing to the disc midplane and r cy \ is the radial distance from 
the z-axis. The quantity po is determined by the mass of 




Fig. 1. Sketch of the general model setup. Shown are the 
scale heights h sma n and /iiarge of the two dust populations 
with h = ^ 3ma " . Note that the scale height is not the upper 



border of a dust population. 

the entire disc, i?ioo = 100 AU is the radial normalisation 
constant, and h the local scale height: 



h(r cyi) = h 



r cy i 



R 



100 



p 



(4) 



Here, a, (3, and ho in Eqn. (j4|) are geometrical parame 
ters. This disc model has served well already in des crib 



ing observed discs, such as HH 30 (IWood et al|l2002l) . the 
Butt erfly-Star IRAS 04302+ 2247 (IWolf et all l2003l). HK 
Tau (IStapelfeldt et al.lll998D. IM Lupi dPinte et al.ll2008h 



HVTau dstapelfeldt et alJl2003ft . and CB 26 dSauter et all 
12009ft . 

To incorporate dust settling and coagulation, this 
ansatz has to be extended. A considerable amount of 
work has been done to understand the physics of dust 
grain coagulation and settling, including extensive nu- 
merical simulations on the topic, as well as dust coag- 
ulation studies in laboratories. A comprehe nsive treat- 
ment c an be found in a s eries of papers by Blum et al.l 
(120061): iLangkowski et all (j2008ft : IWeidling et all (|2009h : 
iGiittler et al.1 (|2009h . Detailed models of dust grain growth 
and settling and the resulting vert i cal di sc structure can 
be found in iDullemond fc Dominikl (|2004t ) . An analysis of 
the ra dial drift of dust gra i ns can be found in lBrauer et al.l 
(|2008ft and iBirnstiel et all (|2010D . However, in the present 
work, we do not include the detailed picture of the 
dust growth process developed there. Instead, we use the 
parametrised ans atz that has been deve loped and success- 
fully employed bv lD'Alessio et all (|2006ft . This enables us to 
formulate a parameter space that includes not only basic 
parameters of the coagulation/sedimentation process but 
also important observational parameters that do not in- 
fluence the disc physics but the way observations perceive 
them. 

To describe grain growth and settling, our model fea- 
tures two distinct dust populations which differ only in the 
maximum grain size. For each population, the minimum 
grain size is the same as found in the ISM: a m ; n = 5 nm. The 
first population has the same maximum grain size as com- 
monly assumed for the interstellar medium: a max = 250 nm 
(jMathis et allll977l) . This dust population is considered the 
reservoir of small dust grains from which the second popula- 
tion of larger grains is fed. The second component features 
a maximum grain size of 1 mm. 
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As described by iDullemond fe Dominikl (|2004l) . the 
height to which a dust grain settles does not only depend 
on its size but also on the level of turbulence present in the 
disc. To model dust settling both grain populations feature 
different scale heights that are related via 

7 h>Q small grains /-\ 

h=— g . (5) 

"'O, large grains 

Here, /io,smaii g rains and Viargograins are the respective scale 
heights. Figure [T] depicts the model concept. 

In the initial state both populations have the same dust- 
to-gas ratio. However, once the disc evolves towards the 
grown and settled case, the small dust grain population gets 
depleted and mass is tr ansferred to the l arge d ust grain pop- 
ulation. The model of ID'Alessio et all (|2006l ) parametrises 
the relative weighting of the two populations by the deple- 
tion parameter e which is the relation of the dust-to-gas 
ratio in the small grain population Csmaii to the default 
dust-to-gas ratio C = iY Gas = 100: 

Csmall / n \ 

e =— — ■ (6) 

As no accretion by the central star nor mass loss to the 
ambient space is considered, the dust-to-gas-ratio Ciargc of 
the large grain population has to be altered in order to keep 
C at its initial value. The values of C iaree a s a function of 
Csmaii are taken from ID'Alessio et"al~l (|2006|) . Also, we keep 
e independent of the radial position in the disc. 

2.2. The parameter space 

Our model features among intrinsic parameters the observ- 
ing wavelength and the disc inclination as external param- 
eters. We now discuss their value and range. 

2.2.1. Stellar heating 

The stellar radiation heats the dust which then itself re- 
emits at longer wavelengths. We neglect accretion and tur- 
bulent processes within the disc as another possible intrin- 
sic energy source. In this sense the disc in our model is 
passive. We e mploy an "avera g e" TT auri star as given by 
the survey of iGullbring et all (|1998t) . This star has a ra- 
dius of r = 2R Q and an effective surface temperature of 
T er j = 4000 K. Assuming the star to be a blackbody radia- 
tor, this yields a luminosity of = 0.92L Q . 

2.2.2. Dust properties 

We assume that the gas within the disc is optically thin 
and limit ourselves to radiative transfer through the dust 
and neglect line emission and absorption by the gas. This 
is justified as the dust is the dominant coolant of the disc. 
Hence, thermal equilibrium obtained in radiative transfer 
calculations based on dust only will provide a reliable de- 
scription of the disc's th ermal structure. This assum ption 
is based on the work of iKamp fe Dullemondl (|2004f ) . who 
showed, that dust and gas temperature in a circumstellar 
disk around a T Tauri star are well coupled in the regions 
which arc sufficiently dense to be traced with current and 
near-future (sub)millimctre observatories. 



We assume the dust grains to be homogeneous spheres. 
Real dust grains are expected to feature a much more com- 
plex a nd fractal structure. As discussed by IVoshchinnikovl 
(|2002f ) , chemical composition, size and shape of dust grains 
cannot be determined separately, but only as a combina- 
tion. 

We keep the dust composition of silicate graphite 
material fixed in order to avoid degeneracies. For the 
optical data we use the complex refractive indices of 
"s moothed astronomical si li cate" and graphite as published 
by IWeingartner fc Draind (|2001h . For graphite we adopt 
th e common — |" appro ximation. As has been shown 
by iDraine fc~M alhotral (|1993|) , this graphite model is suf- 
ficient for extinction curve modelling. Applying an abun- 
dance ratio from silicate to graphite of 1 x 10~ 27 cm 3 H _1 : 
1.69 x 10 _27 cm 3 H~ 1 , we get relative abundances of 62.5% 
for astronomical silicate and 37.5% graphite (^ey and §£_!_)• 
As the grain size distribution we assume a power law of the 
form 

n(a) da ~ a -3 ' 5 da with a m ; n < a < a max . (7) 

Here, a is the dust grain radius and n(a) the number of dust 
grains with a specific radius. For each dust grain ensemble, 
1000 logarithmically cquidistantly distributed grain sizes 
within the interval [a m i n : a max ] have been taken into ac- 
count. For simplicity reasons, we assume Eqn.[7]to be valid 
in the whole dis c for each grain population individually. 
iKim et all (|l994l ) argued that in order to achieve good fit to 
U, B, and V linear polarisation measurements, the grain size 
distribution has to depart from the power law. However, the 
grain size distribution in a specific region of the circumstel- 
lar disc depends on the grain evolution there which depends 
on a variety of factors. First, grains are understood to grow 
fast in the denser regions of the disc. Secondly, mixing and 
fragmentation processes will also influence the grain size 
distribution. With the model employed in this paper util- 
ising two different grain populations we consider effects of 
the vertical dust grain settling. The overall grain size dis- 
tribution then deviates from a simple power law accord- 
ingly. Finally, we assume an average grain mass density of 

Pgrain = 2.5 g Cm" 3 . 

2.2.3. Parametrisation of growth and settling: e and h 

In the framew ork of our model we fo llow the parametric 
values given bv lD'Alessio et all (|2006t) for the depletion pa- 
rameter e and the relative sedimentation height h. For co- 
agulation and settling induced deviations from the default 
mass-to-gas relation of the small grain population we have 
e = 1.0, 0.1, 0.01, and 0.001. The quantity h describes the 
level of turbulence in the disc. As this level is generally not 
known, we adopt three values for ho, scattered around the 
ID'Alessio et~aT1 (|2006l) value of 10. Namely we have h = 8, 
12, and 20. 

2.2.4. Disc flaring and surface density 

In our model (Eqn. [3]) two exponents a and j3 control the 
density distribution of the disc. The exponent p of the 
surface density £ = Y,or p is directly related to a and (3 
via p = a — p. Most observed circumstellar discs can be 
well described by a relatively small range of values for a 
and /3. Good agreement between model and observations 
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are found for values of /3 = | and for a by the relation 
a = 3 (P-k), which is a r esult of viscous accretion theory 
(|Shakura fc Sunvaevll 19731) . Thus, these two parameters are 
also kept fixed in order to reduce the free parameters of the 
mod el. We employ the values as obtained for the Butterfly 
star dWolf et al.l l200l: a = 2.37 and /3 = 1.29. 



2.2.5. Inner and outer radial extent of the circumstellar disc 

Earlier and recent modelling efforts have shown a variety 
of size scales of the inner disc radius. Some discs are best 
described using an inner radius as small as the sublima- 
tion radius for silicate at T — 1500 K. This temperature is 
reached at a distance of about 0.07 AU from the T Tauri 
star in our model. An examp le for this configuration can 
be found in I Wolf et ail (|2003l ). Other discs, such as HH 30 

apparently feature a large inner 
45 AU from the star. To incor- 



(jGuilloteau et al.l 120081 ) 



void with radii up to 
porate both types of discs, we allow the inner radius to 
assume the values ri nnor = 0.1 AU, 1.0 AU, and 10 AU. The 
outer radius of the dust disc is also subject to variation. 
For discs around T Tau ri stars, disc ra dii vary from 25 AU 
(jNeuhauser et al.ll2009t) up to 400 AU ([Simon et al.ll2000[) . 
Hence, we allow in our study the outer disc radius to assume 
the values r out0 r = 100 AU, 300 AU, and 400 AU. 



2.2.6. The total disc dust mass 

The values for the total dust mass of the disc span a range 
of three orders of magnitude: m^st = 3.2 x 10 _6 M Q , 3.2 x 
1O- 5 M , and 3.2 x 1O- 4 M . 



2.2.7. Choices of inclination 

The inclination under which a circumstellar disc is observed 
is a non-intrinsic quantity of the model. In the case of face- 
on discs (8 — 0°) one has direct observational access to the 
central star while in the edge-on case the star is typically 
completely obscured by the dust disc. On the other hand, 
the edge-on case might allow to examine the vertical struc- 
ture better if upper disc layers are sufficiently transparent. 
It is one goal of this study, to determine to which extent the 
vertical disc structure can be determined by observations. 
We include both extrema in our study. 



2.2.8. Observing wavelengths 

To predict observational quantities in a multi-wavelength 
set-up, we compute images in the near-infrared, in the 
mid-infrared, and in the sub-millimetre and millimetre 
regime. Altogether, this gives us images in the I, J, 

H, K, L, M, N, and Q Band (central wavelength A = 

I. /im, 1.25 /im, 1.6 /im, 2.2 /jm, 3.4 /jm, 4.7 /im, A = 10 /im, 
and 20 /im) as well as 15 /im, 50 /im, 100 /im, 350 /im, 
670 /im, 840 /im, 1100 /im, and 1300 /xm. For computations 
of observed fluxes a distance of 140 pc is assumed, repre- 
sentative for the well-studied low-mass star-forming region 
in the Taurus constellation. 



2.3. Radiative transfer simulations 

For the continuum radiative transfer simulat i ons we mak e 
use of the program MC3E0 (|Wolf et al.lfl999l: IWolj 120031) . 
It is based on the Monte-Carlo method and solves the 
continuum radiative transfer problem self-consistently. The 
program u ses the temperature correc tion technique as de- 
scribed bv lBiorkman fc Wood (120011 ). the absorption con- 
cept as introduced bv lLucvl (I1999J) and the enforced scat - 
tering scheme as proposed bv ICashwell &: Everett! (|1959( ). 
Multiple and anisotropic scattering is considered. The ra- 
diative transfer is simulated at 100 wavelengths, logarith- 
mically distributed in the wavelength range [A m j n , A max ] = 
[50 nm, 2.0 mm]. 



2.4. Beyond the model 

Our model setup is based on several simplifications. The 
shape of the disc, governed by a and /3 in Eqs. [3] and 0] is 
kept fixed. Flaring of the disc surface given by /3 has a sig- 
nificant influence on the amount of intercepted stellar radia- 
tion, consequently the thermal structure of the disc, as well 
as the spatial distribut ion of scattered light. As discussed 
by IWatson et al.1 (|2007t) , this two parameters have also an 
impact on the width of the dark lane similar as depletion 
parameter e (see Sec. 13. 1|) . Consequently, the freedom to 
vary a and /3 harbours a degeneracy in this respect. Yet, 
this is not necessarily a draw back: In order to understand 
the influence of all possible physical processes, careful in- 
vestigation of each is needed separately. To do this for the 
effects of sedimentation on images in the described frame- 
work is the goal of this paper. 



3. Results 

We investigate the parameter space of our model for ef- 
fects on various observational quantities that allow one to 
constrain dust settling and coagulation. In particular, we 
examine the implications of the depletion parameter e and 
the scale height relation h in Section 13.11 The general re- 
sult is that dust settling and coagulation as described by 
the depletion parameter e can be constrained by various ob- 
servations if the disc is seen edge-on. The respective find- 
ings are presented in Section 13.21 and 13.31 Possibilities to 
constrain the relative scale height h (i. e. the amount of 
turbulence in the disc), is discussed in Section l3~4l These 
findings hold true only in the parameter space outlined in 
Section 2, especially the exact value of modelled fluxes or 
the morphology of images, particularly in the mid-infrared. 
Nevertheless, the qualitative behaviour of tracers for grain 
growth and sedimentation outlined in this Section can be 
also expected for models with, for instance, different choices 
for a, j3 or ho in the density distribution (cf. Eqn. [3]). 

3.1. The fiducial model 

To discuss the effects of dust grain growth and coagulation 
on observable quantities and as a reference for further inves- 
tigations we introduce a fiducial model. Its parameters are: 
rw = 0.1 AU, r outor = 300 AU, m dust = 3.2 x 10" 4 M o , 
h = 8, and e = 1.0. This approach allows us to evaluate 



1 Available upon request: wolf@astrophysik.uni-kiel.de 
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Fig. 2. Dependency of the width of the dark lane in scat- 
tered light images on the depletion parameter e (grey) and 
the observing wavelength A (black). Shown are cuts through 
images along the disc midplanc centred on the projected 
position of the embedded star. 



the impact of variations of the individual model param- 
eters on the observable quantities. Resulting degeneracies 
are discussed. 

Effects of the depletion parameter e The depletion parame- 
ter e governs how much mass is transferred from the small 
to the large dust grain population. Therefore, e can be un- 
derstood as a parameter for dust grain coagulation and the 
resulting sedimentation. 

The near-infrared images of edgc-on discs exhibit a pro- 
nounced chromaticity. In this wavelength regime the disc's 
midplane is optically thick and shields the central star from 
direct observation. The height above the midplane at which 
the disc becomes transparent to the scattered stellar light 
constrains the disc's scale height. 

Besides the optical properties of the dust particles, the 
opacity of the disc depends on the amount of dust from 
both dust populations present in the respective disc layers. 
Considering dust settling, the composition of those layers is 
a function of the degree of depletion. Thus, the width of the 
dark lane is not a function of the disc gas scale height alone 
but must also be a function the depletion parameter e. The 
fiducial model exhibits this behaviour. Fig. [5] shows both, 
the dependency of the dark lane on the observing wave- 
length A as well as on the depletion parameter e. Although 
the values for A and e in the plot do not exactly mirror 
each other the width of the dark lane can be reproduced 
adjusting both parameters. 

Fig. [3] shows the relation of the radial brightness profile 
of the fiducial model with different total dust masses and 
different depletion parameters e observed face-on. For r = 
0, the ratio between the fluxes of the fiducial model with e = 
1.0 and with e = 0.001 1 since the same stellar properties 
are used for all models. For r > 0, the ratio between the 
brightness profiles is « 2 for different depletion parameters 
e while it is slightly sloped for different total dust masses. 
A variation of the density structure (i.e. a and (3 in Eqn. 
[3]) can easily counterbalance the effects shown in the plot. 
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Fig. 3. Selected flux ratios of the radial brightness profile 
of the fiducial model seen in the H band face-on. The solid 
line shows the flux ratio between models with e = 1.0 and 
e = 0.001. The dashed line shows the ratio between fiducial 
models with mdust = 3.2 x 1O _4 M0 and mdust = 3.2 x 
10- 6 M ra . 
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Fig. 4. Optical depth of the fiducial model along radial 
paths for selected inclinations as a function of wavelength. 
A high optical depth in the millimetre regime is obtained 
only for high inclinations while for shorter wavelengths the 
model is opaque in a wider range. 

For inclinations 9 = 60° and 9 = 0° the dependence of the 
brightness distribution on e is almost lost. 

The first reason for this ambiguity is that as the disc is 
tilted to smaller inclinations, the dark lane vanishes. This is 
because the optical depth along the line of sight decreases 
and more scattered stellar light reaches the observer. 

Secondly, as the disc itself is still optically thick in the 
near infrared, it is still only the upper disc layers that are 
penetrated by stellar light. Hence, it is the same physical 
environment as in the edge-on case that determines the 
characteristics of the scattered light. 

Fig. [4] shows the dependency of the optical depth as 
a function of the wavelength for selected radial paths. A 
high optical depth (t\ 3> 1) in the millimetre regime is 
obtained only at high inclinations. For shorter wavelengths, 
especially in the ncar-infrarcd, the fiducial model is opaque 
even for 9 ~ 70°. Figures [5] and H2 show the dependency of 



J. Sauter and S. Wolf: Dust grain growth and settling in circumstellar discs 



7 



10" . 




10000 



Fig. 5. Optical depth of the fiducial model along a radial 
path at 90° inclination as a function of wavelength. Only a 
weak dependency on the depletion parameter e exists. 
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Fig. 6. Optical depth of the fiducial model along a radial 
path at 85° inclination as a function of wavelength. Unlike 
the situation in the exact edge-on case (Fig. [5]), the optical 
depth becomes sensitive to the depletion parameter e by 
typically half an order of magnitude. 



T\ along a radial path on the depletion parameter e. For the 
edge-on orientation 6 = 90° the dependency t\ (c) is very 
weak (Fig. [5]). Yet, 5° away from the edge-on orientation, 
the optical depth is much more sensitive to the depletion 
parameter e. 

For images taken in the sub-millimetre and millimetre 
regime no distinctive effect is observed. While absolute flux 
values change as the large dust grain population becomes 
more massive, no relative effect is seen either with respect 
to the location within one image or as a function of wave- 
length. 

These findings hold true for any inclination we consid- 
ered in our framework. Since for smaller values of e more 
large grains are present in the system and large grains re- 
emit more efficiently at larger wavelengths than smaller 
grains do, one might expect observable effects in the mil- 
limetre regime. For the case of an edge-on disc the results 
are discussed in Section I3.2I However, if the disc is seen 
face-on, there are no observable effects. 
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Fig. 7. Contour lines of the density distribution normalised 
to the respective peak density of the two grain populations. 
The two panels show the effect of small (top) and large 
(bottom) relative scale height h. 

Effects of the height relation h The height relation h as de- 
fined in Eqn. [3] describes the level of turbulence in the disc 
as well as the magnitude of the viscosity parameter a. Fig. [7] 
illustrates the effects of h on the density distribution. The 
density distribution of the small grain population does not 
depend on h. According to the model set-up the large grain 
population becomes concentrated in a much smaller space 
for larger h. 

The effects of h on images at different inclinations is 
essentially that of an amplifier of effects caused by the de- 
pletion parameter e. In the case of near-infrared images of 
edge-on discs, the width of the dark, high-opacity lane is 
affected. Larger values for the height relation h result in a 
dark lane with larger contrast and vice versa. This can be 
readily understood in terms of optical depth. Higher values 
for h essentially squeezes the large dust grain population 
into a smaller volume and thus increases the optical depth. 
The amplifying effect is stronger for smaller values of e: 
more mass is packed into the large dust grain population 
yielding a height optical depth. The magnitude of the ef- 
fects caused by the variation of h is considerably small. 

In Fig. \8\ the optical depth of the fiducial model along 
a radial path in the edge-on orientation as a function of 
wavelength is shown. The comparison with Figs. [5] and |B] 
shows that for h the dependency is not as strong as for e. 

3.2. Millimetre Images 

Millimetre images are well suited to classify dust settling 
and coagulation because at these wavelengths the dominant 
fraction of the radiation has its origin in the thermal re- 
emission of the dust grains. Further, the redistribution of 
the mass also alters the optical depth structure. 

The most illustrative way to see how millimetre images 
provide means to constrain disc evolution (the changing ra- 
dial and vertical opacity and thus the reemission structure) 
is to compare different regions of a disc image, especially if 
the disc is seen edge-on. The regions of our choice are de- 
picted in Fig. I^and labelled (cw, starting in the image/disc 
centre) A, B, C, and D. One expects that as the disc evolves, 
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Fig. 8. Optical depth of the fiducial model along a radial 
path at 90° inclination as a function of wavelength. Plotted 
is the dependence on the relative scale height h. A compar- 
ison with Fig. |5] shows, that the dependence on h is not as 
pronounced as the dependency on e. 
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Fig. 9. Location of the four regions used to compare in- 
tensities. The dashed lines are logarithmic contour lines of 
the brightness distribution of the fiducial model observed 
at A = 1.3mm and 9 = 60°. 



the layers distant from the midplane loose more and more 
mass due to the dust grain growth and settling towards the 
midplane, the amount of thermal re-emission from those 
layers decreases as well. Radiation originating close to the 
midplane is expected to raise or to stay at the same level, 
since the growing midplane mass not only yields more emit- 
ting grains but also more absorbing grains (important in 
the submillimcter regime) . Densities values are already or- 
ders of magnitude higher close to the disc midplane, the 
down-raining grains therefore cannot significantly increase 
the optical depth and thus level any gain in emitter num- 
bers. 

Fig. [TU] shows excmplarily our findings by relating the 
flux from the four different regions to the total flux. Starting 
from red for the early, e = 1.0 case and ending at blue with 
a depletion parameter e = 0.001. The shape of any mark de- 
notes the region from whence it is taken. Plotted is the rela- 
tive flux from a certain region versus the size of that region. 
The latter can be understood as the degree of resolution of 
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Beamsize [3 AU] 
Wavelength= 1.3 mm, Inclinations 90* 

Fig. 10. Analysis of image brightness distribution of a cir- 
cumstellar disk at different evolutionary stages. The results 
for different evolutionary states are coded in colour. The 
shape of any mark denotes the region from whence it is 
taken. Vertical lines indicate resolution limits of selected 
interferometer. Horizontal lines give the ALMA sensitivity 
limit, the other two instruments do not have the required 
sensitivity. The flux is given in units of the total flux. Since 
this flux depends on e as well, the sensitivity limit is af- 
fected. The dependency of the flux in regions B and C on 
e is shown. 



the observing instruments. As an indicator, the maximum 
resolution 1.3 mm of various instruments (ALMA, PdBI, 
SMA) is given as dashed black lines. Although all three in- 
terferometers shown in the plot provide the required resolu- 
tion to detect the effect, only ALMA provides the required 
sensitivity. As the flux is given in units of the total flux 
of the image, which depends on e, the sensitivity is also a 
function of e. As expected, fluxes taken in the midplane do 
not vary considerably with the evolutionary state: all four 
colours are almost atop of each other. However, for fluxes 
taken from above the disc midplane, a clear decrease in flux 
is visible. From the initial e = 1.0 to the final e = 0.001, it 
spans almost two orders of magnitude. This effect of flux 
decrease though is limited to the edge-on case, only. 

At inclinations of 9 = 60°, 0°, the lines of sight are no 
longer within the midplane or any parallel plane but inter- 
sect higher disc layers as well as the disc midplane. As a 
result the radiation comes from various disc layers, repre- 
senting different stages of dust coagulation and sedimen- 
tation. The amount of flux collected in the four different 
regions is not a function of the depletion anymore but just 
of the disc's geometrical thickness. As the inclination can 
be determined from image-morphology, these findings are a 
valuable tracer for dust sedimentation in the edge-on case. 
To estimate the robustness of this result, we consider the 
influence of disc parameters other than the depletion factor 
e. We find that although the total flux value from different 
regions changes, the effect of separation of the flux values 
in the upper and lower layers of the image with later evo- 
lutionary states does not depend on the chosen parameter 
value. It is also independent of the observing wavelength in 
the (sub) millimetre regime. 
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Fig. 11. Different morphologies of an edge-on disc in the 
N band with r ilmcr = 0.1 AU. Top left: m dust = 3.2 x 10~ 6 , 
Pouter = 300 AU, h = 12, e = 0.001. Top right: mdust = 
3.2 x 10~ 4 , r out0 r = 300 AU, h = 20, e = 0.1. Bottom left: 
m du8 t = 3.2 x 10~ 4 , r outC r = 300 AU, h = 20, e = 0.1. 
Bottom right: mdust = 3.2 x 10 -4 , r ou ter = 400 AU, h = 8, 
e = 1.0. All images feature a linear colour scale. 

3.3. The disc in N and Q band edge-on 

We now turn our attention towards the potential of mid- 
infrared observations in the N and Q bands. In our cho- 
sen models, the re-emission in these bands dominates 
the scattered stellar flux by several orders of magnitude. 
Morphologically, it is the transition region from the 'dark 
lane' type of brightness distribution as seen in the near- 
infrared (see Section 14. 2\i to the extended and elongated 
typical for (sub) millimetre images. In our edge-on images 
the typical central dark lane is very thin, in the order of 
one resolution element of the Giant Magellan Telescope. 
The major difference to millimetre images of edge-on discs 
is that the flux due to re-emission is concentrated in the cen- 
tral region while at increasing wavelengths the emission be- 
comes increasingly extended as cooler material dominates 
the emission. 

Based on our model we predict various morphologies 
that can be observed in this wavelength regime. Possible 
scenarios from within our parameter range include a sin- 
gle central peak, two elongated emission regions separated 
by a dark lane, three emission maxima in the midplane, a 
four-legged x-like shape and others more. Fig. [TT] gives an 
impression of the variety of possibilities. The reason for this 
zoo of different morphologies resides in the fine-tuning of 
the optical depth. Depending on the total mass and its dis- 
tribution, different regions of the disc are rendered opaque 
enough for a specific choice of values parameters such that 
their radiation contributes to an image (or not). In most 
cases t\ is close to unity within an order of magnitude. 

However, the overall width of the brightness distribu- 
tion turns out to be a valuable tracer for dust grain growth 
and sedimentation in the edge-on case (see Fig. I12p. A lin- 
ear relation between the FWHM of the emissions in the 
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Fig. 12. FWHM of the height of the emission in N band 
versus the same in Q band for different values of the de- 
pletion parameter e. A general tendency of a linear relation 
between the FWHMs for smaller values of e can be ob- 
served. 



Depletion e 


linear slope 




1.0 


1.0 


10715 


0.1 


0.9 


1820 


0.01 


0.8 


381 


0.001 


0.8 


191 



Table 1. Linear slope and % 2 of the N/Q band width of 
the models. 



two bands is observed. For all possible combinations of pa- 
rameters in our framework, the scattering of the Q band 
peak width vs. the N band peak width is strongest for a 
large depletion parameter and ceases for later evolutionary 
stages. 

Introducing a straight line as a fit, the chi-square (x 2 ) 
is a function of the depletion parameter e as shown in 
Table Q] The reason for this increasing linearity is the 
fine-tuning of the optical depth. Re-emitted radiation re- 
ceived from the object originates at regions with an optical 
depth smaller than unity. For large values of e the disc 
dust mass is distributed over a larger volume. Thus, the 
iso-surface with TNband has more morphological freedom. 
Consequently, the width of the brightness distribution of 
one disc model picked at random depends more on the set- 
ting of the other parameters than this is the case for smaller 
values of e. There the bulk of matter is concentrated close 
to the midplane and the shape of the regions seen is less 
subject to a specific choice of parameters: The density dis- 
tribution of the large grain population determines the width 
of the brightness distribution alone, not the interplay be- 
tween two grain populations. The specific value of the slope 
in Table [1] is roughly determined by the wavelength depen- 
dency of r and the vertical density increase of the large 
dust grain population as given by Eqn. [21 Increasing the 
observing wavelength from 10 (im to 20 fim allows a deeper 
look into regions of higher disc densities. For e = 0.001 a 
slope of 0.8 implies that the Q band image is a little thinner 
than the N band image. This means that the spatial gain 
by reduced optical depth is almost counterbalanced by the 
increased density in deeper layers. 
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Within our parameter space the FWHM-relation of the 
emission thickness in N and Q band for one disc depends on 
the total dust mass. An analysis of disc models shows, that 
they group in Fig. [T2] according to their total dust mass, 
not their relative scale height. This is also true for later 
evolutionary stages. 

Consequently, the evolutionary stage of a single ob- 
served disc can be inferred by comparing the extent of the 
emission in the N and Q band only. The precise value of 
this relation will depend on the global geometry of the disc. 
Nevertheless, a peak width ratio of ~ 0.8 from Q to N band 
is a good indicator for discs with most of their masses set- 
tled in large dust grains close to the midplane. 



3.4. Q and N band face-on: Turbulence digs in 

This Section is dedicated to investigate the possibility to 
constrain the relative scale height h. Fig. [TT] already indi- 
cates that images in the N and Q band are sensitive to the 
fine-tuning of the density distribution. Fig. [T3] illustrates 
this in the N band for the case of face-on discs. In this 
regime the observed flux originates predominantly from re- 
emission of hot dust at ~ 300 K. The marker for the relative 
scale height is the gap in the flux distribution at intermedi- 
ate disc radii of ~ 100 AU. The strength of this emission gap 
depends directly on the relative scale height h. For large h, 
i. e. a small scale height of the large dust grain population, 
the gap is clearly visible. It ceases for smaller values of h. 
The effect also depends on the observing wavelength. 

The nature of this apparent gap can be understood 
by analysing the spatial temperature distribution together 
with the optical depth. Fig. 1 141 shows temperature contour 
lines of the large-grain population enclosing regions with 
temperatures higher than the black-body emission peak in 
the Q and N band. This emission peak occurs for the Q 
band at 144 K and in the N band at 315 K. Regions out- 
side are cooler. Closer to the midplane the reason for this 
is the high optical depth such that the material there can- 
not be heated effectively by stellar radiation. Far above 
the midplane the density of the disc is too low to assign 
a meaningful temperature (cf. Fig. [7]). Fig. RBI also shows 
the TQ t N = 1 isoline for the two bands obtained by inte- 
grating perpendicular to the midplane, starting at z = oo. 
For both bands and the two extreme values for h, the high 
temperature regions of the small grain population are far 
above the tq n = 1 line and are not shown in the plot. In 
both bands, these regions most efficiently contribute to the 
face-on image of the disc for any value of h. These regions 
are the origin of the flux peak in the centre, cf. Fig. [13] 

In the Q band also the large-grain population con- 
tributes to the image for both values of h. Unconvolved face- 
on images show that the surface brightness of the small- 
grain population is about an order of magnitude larger. 
The emission region of the small grain population is close 
to the rotation axis of the system. Seen from above, all 
emission of the small grain population is seen piling up in a 
narrow region while the emission regions of the large grain 
population are stretched above the midplane. In the edge- 
on case, this geometrical effect yields the opposite result: 
The emission from the large grain population is seen in a 
smaller area and thus contributes more to the image than 
the small grain population. 
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Fig. 13. Face-on images in the N band and radial profiles 
of the brightness distribution. Top: Fiducial model with a 
scale height relation of h = 8. Centre: Fiducial model with 
a scale height relation of h = 20. Bottom: Radial profiles. 
The solid lines corresponds to the first image, the dashed 
line to the second. All images have been convolved with a 
PSF of 30 AU. The colour scale is chosen as - S?; 1 . 



Matters arc different in the N band. An optical depth 
of tn = 1 is reached at higher altitudes from the mid- 
plane than in the Q band. For h = 20 and small radii, 
the emitting region is below the tq = 1 isoline. This leads 
to a more dimmed emission in the N band for middle radii 
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Fig. 14. Temperature distribution and isolines of constant 
optical depth r = 1 in the N band (top) and Q band (bot- 
tom) for the fiducial model. Different colours refer to differ- 
ent relative scale heights h. Note the inclusion/exclusion of 
the 315K - emission region in the N band of the large-grain 
population 
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Fig. 15. Classification of the apparent gap in the bright- 
ness distribution of edge-on seen in the N band. Shown is 
the total brightness minimum in units of the average disc 
brightness in the outer 73 of the disc versus the minimum 
brightness in the outer 73 of the disc. Colour coded are 
model sets with different parameter values. 



which results in the brightness gap at these radii. The iso- 
lines are just indicating how deep face-on observations pen- 
etrate into the disc. Wien's law only allows to determine 
those grains and disc layers which emit most effici ently in 
a given wavelength band (e.g. IScheeerer et al. 2006]). Many 
more grains with lower or higher temperature can also pro- 
vide the dominating flux at the wavelengths in question. 
However, in Fig. all regions with higher temperatures 
are included by the temperature contour lines. Considering 
that the emitted energy is proportional to ~ T 4 , the rea- 
soning in the previous paragraph remains robust. 

Fig. [TS] shows that the contrast of the brightness gap is 
a function of the relative scale height for two thirds of our 
parameter space. Plotted is the total brightness minimum 
in units of the average disc brightness in the outer third of 
the disc versus the minimum brightness in the outer third 
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Fig. 16. Classification of the brightness gap of face-on 
seen discs at A = 10 yum. Plotted on a log-linear scale is 
the total brightness minimum in units of the average disc 
brightness in the outer / 3 of the disc versus the gap width. 
Colour coded are models with different scale height relation 
h. Models featuring h = 8 do not exhibit a gap, see Fig.[T5l 

of the disc also in units of the average outer disc brightness. 
Models featuring a gap have a flux minimum within the first 
two thirds of the disc and are located left of the line with 
unity slope. Models without an inner gap appear on the line 
with unity slope. The gap is typically located at radii r rj 
100 AU. Discs with an outer disc radius of r ou t er = 100 AU 
consequently do not exhibit a gap. In Fig.[15l these discs are 
all located on the line with unity slope. We define the width 
of the gap as the radial distance between the location of 
the first brightness maximum after the gap and the location 
where the brightness distribution close to the centre reaches 
the maximum value. Fig. [If)] shows the minimum flux in 
the gap as a function of this gap width. Models with very 
pronounced gaps are thus located in the lower right part of 
the plot. Shown are models with different values of h. The 
strength of the gap does not correlate with h as can be seen 
from Fig. [T5] Yet, the width of the gap docs. Our analysis 
further shows that the gap width does not correlate with 
the outer disc radius r out0 r. 

4. Other markers: Discussion 

Our study also shows, that markers others than those pre- 
sented in Section [3] are not sufficient to constrain grain 
growth and coagulation. However, we discuss them in this 
section to shed further light on the nature of the images of 
circumstellar discs. 

4.1. Radial profiles in the millimetre regime at low 
inclinations 

The observed flux in (sub-)millimetre images stems from 
dust grain re-emission at these long wavelengths. The op- 
tical depth is small rendering all dust particles visible. As 
our models are rotationally symmetric, it suffices to con- 
sider the radial brightness profile instead of the whole im- 
age. However, those profiles do not provide means by which 
discs at different stages of their evolution could be distin- 
guished. While the absolute value of the radial profile of the 
surface brightness distribution depends on the evolutionary 
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Fig. 17. The chromaticity and dependency of the width of 
the dark lane on the depletion of upper disk layers in the 
NIR wavelength regime. For longer wavelengths and smaller 
e the width is generally smaller. Shown also are straight line 
fits. 



stage of the disc, the shape of the profile does not. This ren- 
ders the radial profile for 8 < 90° degenerate with respect 
to TOdust and e. This behaviour agrees with the findings of 
Section [~~ 



4.2. Dust lane chromaticity in the near-infrared 

In the near-infrared, the typical appearance of a circum- 
stellar disc seen edge-on is a bipolar structure that is inter- 
sected by a dark lane. The width of the dark lane is an effect 
of optical depth and wavelength dependent. In addition, 
at wavelengths between 2 /im and 10 fim, radiation due to 
thermal re-emission of the dust grains becomes increasingly 
important in our disc models and is in the mid-infrared by 
far the dominant contributor to the total flux. As the op- 
tical depth is a function of the disc mass, one expects the 
width of the dark lane to be dependent also on the deple- 
tion of the upper disc layers. As grains grow and rain down 
towards the midplane, the upper layers become less opaque 
and might affect the width of the dark lane as an indicator 
for e. Increasing the observing wavelength unveils deeper 
layers in the disc and the large dust gain population (cf. 
Fig. W) increasingly influences the appearance of the disk. 

Fig. H7J shows the width of the dark lane for the fiducial 
model. As expected, the width of the dark lane decreases 
with both, wavelength and with progressing depletion of the 
upper disc layers. The chromaticity itself is a function of e 
(otherwise the fitted lines would be parallel): For smaller e 
the chromaticity of the dark lane is less steep. 

Moreover, the chromaticity's dependence on e provides 
means to identify dust grain growth and settling uniquely 
via near- to mid-infrared imaging only. In Fig. [18] all chro- 
maticity slopes m e for all possible parameter combinations 
within our framework are given. For any given e, the slope 
obtained is either zero or scatters within a certain range. 
The depletion of the upper disc layers could be inferred 
from the slope if a unique relation between e and the chro- 
maticity slope could be observed. This is not the case. The 
same behaviour is observed when not only images in the 
I, J, H, and K bands but also in the L and M bands are 



h = 8 
h-20 
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Fig. 18. Scatter plot of the chromaticity-slope of the dark 
lane for different e and different fit ranges. 
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Fig. 19. Width of the dark lane as function of the FWHM 
of the millimetre emission. Colour coded is the depletion 
parameter e. A clear dependency between e and the width 
of both is not seen. 



considered. These two bands mark the transition region be- 
tween images dominated by scattered stellar light and dust 
grain rcemission. This is independent of the resolution used 
to image a given circumstellar disc in the NIR. 



4.3. Near-infrared and millimetre 

One might expect that a combination of millimetre im- 
ages and near-infrared images yields tracers for dust grain 
growth and sedimentation just as the combination of im- 
ages in the N and Q. Fig. [T2] gives the relation between 
the FWHM of the vertical extend of the spatial brightness 
distribution in the millimetre regime and the width of the 
dark lane in the J Band. The high dust densities in the disc 
mid-plane are responsible for the dark lane in the NIR and 
the elongated spatial brightness distribution in the millime- 
tre regime. Thus, a relation between the width of the dark 
lane and the millimetre emission region could hint to the 
disc's evolution. However, as Fig. Q15] shows, no correlation 
between the two quantities is observable. 
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Fig. 20. The FWHM of the spatial intensity distribution 
of disc observed face-on as function of the wavelength. The 
depiction parameter e is encoded with different symbols. 



4.4. Degeneracy between parameters and panchromatic 
studies 

We point out that if all disc parameters in our parameter- 
space arc known, in most cases the chromaticity of the dark 
lane uniquely determines the evolutionary state, as illus- 
trated by the fiducial model (cf. Section 153)1 . In general, it 
is not possible, though, to determine the evolutionary stage 
of a disc when other parameter values remain unconstrained 
as well. 

As a simple counter-example, consider a system that is 
optically thin at all near-infrared wavelengths, and thus a 
dark lane is not visible and no chromaticity can be seen. 
In our framework, his happens for discs with either a low 
total dust mass, a low scale height relation h, a high deple- 
tion parameter e, or combinations thereof. The total dust 
mass can be constrained by the systems' total flux in the 
millimetre regime where the disc is generally optically thin. 
However, to constrain h a certain total dust mass is re- 
quired. Otherwise the effects of h arc rendered indistin- 
guishable. 

4.5. The general face-on case 

As shown for millimetre images in Section 13.21 the ability 
to constrain dust sedimentation by single imaging is lost for 
6 < 90°. For images in the near- infrared and mid- infrared 
the situation is similar. In Fig.^plots of the FWHM of the 
disc's brightness distribution seen edge-on at different wave- 
lengths arc given. At short wavelengths the emission is very 
small as in the NIR the star is the dominant contributor to 
the image flux. Images exhibit an apparently broader flux 
distribution at larger wavelengths since the origin of the 
observed flux shifts from scattering to re-emission. The size 
of the re-cmission regions growths as the observing wave- 
lengths gets larger and cooler parts become more important 
contributors. Whereas in Fig.[T2]a grouping of data points 
with smaller e to a linear relation between the N and Q band 
emission width can be observed no such relation is seen in 
a NIR/millimetre comparison. Even if a parameter is well 
constrained by some observation no specific behaviour of 



data points with the same depletion parameter e can be 
seen. 

This observation d o es no t contradict the findings of 
iDullemond fc Dominikl ()2004f ). In their work, dust grain 
growth and sedimentation have not been modelled with a 
parametrised model but including radius dependent sedi- 
mentation. In our model framework we do not allow for 
radial depended dust settling. Quicker dust grain growth 
at smaller radii renders the upper disc layers in the in- 
ner disc more transparent and allows to heat the outer 
disc more efficient. The radial brightness profile would be 
shifted and/or spread outward. In contrast, the upper lay- 
ers deplete all at the same rate in our framework. Hence 
the radial brightness profile is unaltered. The latter is only 
true for the relative profile, that is a profile normalised to 
some other disc inherent quantity. For most practical rea- 
sons this would be the total disc flux. The absolute value of 
fluxes changes with a the evolution of the disc (cf. Section 

mj. 

5. Conclusions 

In this study wc explore the effects of dust grain growth and 
sedimentation on images of circumstellar discs at various 
wavelengths and inclinations. As the prediction of observ- 
able quantities of circumstellar discs is highly non-linear 
and harbours several degenerate parameters, we employ a 
simplified, parametrised model. This allows us to discuss 
the general effects of the depletion of upper disc layers by 
grain growth and dust sedimentation using radiative trans- 
fer simulations in the near-infrared, the mid-infrared as well 
as in the (sub-)millimetre regime. In the advent of tele- 
scopes providing high spatial resolution we focus on images 
at these wavelengths thus extending the existing studies on 
SEDs. Although our approach is based on certain simplifi- 
cations, we identify tracers that uniquely allow to constrain 
the evolutionary stage of a disc. 

In particular, our findings have direct implications on 
future observations: 

— The width of the dark lane in NIR images of edge-on 
discs is degenerate with respect the degree of depletion 
of upper disc layers from large grains and observing 
wavelength. Although this wavelength and the degree 
of depletion do not exactly mirror each others impact 
on the width of the dark lane, a given width lane can 
be reproduced adjusting both parameters (see Section 

EH). 

— High spatial resolution also needs to be accompanied by 
high sensitivity in the case of millimetre images. Only 
ALMA will provide the required resolution and sensitiv- 
ity to identify grain sedimentation in millimetre images 
of edge-on seen discs. Settled disc show a relative sur- 
face brightness decrease of ~ 1CP 4 from the mid-plane 
to 150 AU above the disc while non-settled discs show 
a decrease of only ~ 10" 2 (see Section l3T2j) . 

— If the disc is not seen edge-on, then this distinc- 
tion between settled discs and non-settled discs cannot 
achieved by millimetre observations, if the settling is 
independent of the radial position in the disc. 

— A comparison of the FWHM of the brightness distribu- 
tion in N and Q band images allows to constrain the 
degree of dust sedimentation in a circumstellar disc if 
seen edge-on. Evolved discs in our parameter space show 
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a ratio of the Q band FWHM to the N band FWHM 
of ~ 0.8 which depends only very weakly on other disc 
properties. Young discs however show a ratio that de- 
pends strongly on other disc parameters and scatters 
for various combinations of parameter values around the 
average FWHM ratio of ~ 1.0 (see Section [53)1 . 
— The relative sedimentation height h to which parti- 
cles settle can be determined using face-on observa- 
tions in the N and Q band. Due to the complex in- 
teraction between the spatial structure of the optical 
depth and the temperature distribution discs with set- 
tling of large grains close to the mid plane exhibit a gap 
in the observed radial brightness distribution at typi- 
cally ~ 100 AU. However, the effect is only seen in discs 
considerably larger than ~ 100 AU (see Section l3T4|) . 

Furthermore, we discuss ambiguities that prevent to 
achieve the same results for images in other wavelength 
regimes, namely the near-infrared. As a next step a disc 
parametrisation should be chosen that allows to incorporate 
radially dependent sedimentation. Also, modelling of im- 
ages of an individual object will provide a better insight in 
its particular evolutionary stage. The knowledge gained the 
framework of the present paper forms a non-comprehensive 
basis for further analysis of dust grain growth and sedimen- 
tation effects on images. 
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